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Abstract 

Validity of modified finite-size scaling above the upper critical dimension is demonstrated for 
the quantum phase transition whose dynamical critical exponent is z = 2. We consider the N- 
component Bose-Hubbard model, which is exactly solvable and exhibits mean-field type critical 
phenomena in the large- iV limit. The modified finite-size scaling holds exactly in that limit. 
However, the usual procedure, taking the large system-size limit with fixed temperature, does 
not lead to the expected (and correct) mean-field critical behavior due to the limited range of 
applicability of the finite-size scaling form. By quantum Monte Carlo simulation, it is shown that 
the same holds in the case of N = 1. 
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I. INTRODUCTION 



Since the quantum phase transition to Mott insulator from superfluid was observed in 
;he optical lattice system this quantum critical phenomena has been one of hot topics. 
2j This system is effectively described by Bose-Hubbard (BH) Hamiltonian. [3j The zero- 

nnn 

temperature phase diagram of BH model has been well investigated [4J, [5|, laj • There are phase 
transition points called multicritical points whose dynamical critical exponent is z — 1 and 
line of the other type of phase transition called generic transition whose dynamical critical 
exponent is z = 2 on the zero-temperature phase diagram. In this paper, we consider the 
generic transition (i.e., z = 2) in three-dimensional systems. The three dimension (d = 3) 
is above the upper critical dimension d u — 2. Therefore, this phase transition is exactly 
classified and its critical exponents should be identical to those of the mean-field theory. 
To estimate the locations of critical points quantitatively, we frequently apply the finite-size 
scaling to the data of finite-size systems calculated using quantum Monte Carlo (QMC) 
method. 

Above the upper critical dimension, the finite-size scaling (FSS) should be modified due 



to a dangerous irrelevant variable, 
modified finite-size scaling (MFSS) 



7] In contrast to the conventional FSS below d u , the 
8( is not justified by renormalization group or scaling 



theories. However, its validity has been demonstrated for the five dimensional Ising model 
0,0 , [hJ, 0(n) model and 4 model in large- N limit 0, [ijj]. For the quantum phase 
transition with z = 1, below the upper critical dimension, a simple application of the FSS is 
trivially possible by identifying the inverse temperature f3 as just an additional dimension. 
Actually, to estimate the multicritical point quantitatively, Smakov and S0rensen |13j applied 
the FSS with the additional argument (3/L to the multicritical point in d — 2 case where 
the system is below the upper critical dimension because d + z < 4. For the quantum 
phase transition with z / 1, below the upper critical dimension, the application of FSS 
is also possible with the additional argument f3/L z instead of f3/L on the ground that the 
ratio between the correlation time £ T and the correlation length £ to the z-th power is 
£ T /£ 2 = 0(1). [4] Zhao et al, applied the FSS to the case z = 2 and d = 2, which is 
just the upper critical dimension, and succeeded in estimating the phase boundary on the 
zero-temperature phase diagram of their model. [14J] The purpose of the present paper is to 
demonstrate the validity of the MFSS in the case where d > d u and z ^ 1, both by Monte 
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Carlo simulation and by exact solutions. We consider the case z = 2, d = 3, i.e., above the 
upper critical dimension. It seems a natural extension to add the argument 13 1 L 2 to the 

n 

scaling function of MFSS. [15( Namely, we assume that the singular part of the free energy 
F s has the scaling form, 

F s (r, T), P, L) ~ Yp {5L^I\ V L^\ (3/L 2 ) , (1) 

with a universal scaling function Yp, where the definition of the free energy is F = — InH 
with the partition function H, r indicates the coefficient of the term including square of the 
order parameter in the Hamiltonian (e.g., the chemical potential \x or the hopping amplitude 
t in the model ([2]) described below), for 5 indicates the difference from the quantum critical 
point (e.g., 5 = r — r c ), and r] is the field inducing the order parameter. 

The critical exponents for the finite temperature behavior at quantum critical point 
should be identical to those of mean- field theory, e.g., % ~ T~ 3//2 where % is susceptibility. 
However, as shown in Sec. IHIj. the exponents derived by the limit L — > oo of scaling form 
(e.g., x ~ T -5 / 4 ) are different from those of the mean-field theory. The reason of this 
apparent contradiction is that the scaling form (pQ) is valid only when f3/L 2 = 0(1). That 
is, we cannot infinitize L in Eq. ([!]) while keeping (3 finite. In this paper, we show that 
the application of MFSS to the z = 2 quantum critical point is reliable, if the condition of 
validity is satisfied, just as well as the conventional FSS below the upper critical dimension. 

In Sec. [Til we define iV-component BH model. In Sec. IHIt we focus on the N — 1 case 
and show the application of the MFSS to the numerical result of the QMC simulation. In 
Sec. IIVt we focus on the N = oo case, which is exactly solvable even for finite systems, to 
show that the susceptibility obeys the MFSS form under the condition j3/L 2 = 0(1). In 
Sec. El we give a discussion and summary of this paper. 

II. iV-COMPONENT BOSE-HUBBARD MODEL 

We consider the iV-component BH model on the hypercubic lattice whose Hamiltonian 
is described as 

N N j- j- N N 

^ = -^EE + b ^l) - p E E h « h « + 2N £ E E b iAhibai, (2) 

o=l a=l i a=l f}=l i 
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where b^ ai (b ai ) creates (annihilates) a a-type boson at site i, and runs over all pairs 
of nearest-neighbor sites. The symbols t, U, and /x, denote the hopping amplitude, the on- 
site interaction between bosons, and the chemical potential, respectively. The coordination 
number in the hypercubic lattice is Z = 2d. We take the lattice spacing as our unit of 
distance. For concreteness, we consider only three-dimensional case in this paper, (i.e., 
Z — 6.) Generalization to arbitrary dimensions should be straightforward. 
Here, we define the free energy F v as 

F v ee -ilnTr[e-^-" s )], (3) 

TV 
a=l i 

with the field rj inducing the order parameter. 

The iV-component BH model ([2]) is solvable in the large- N limit. In Sec. IIV| we demon- 
strate that the MFSS scaling (CQ) exactly describes the asymptotic behavior of the model 
([2]) in the large- N limit. We note here that an exactly solvable model similar to the present 



one was investigated in the 1980s. [16|, [l7j The model was defined with Bose field opera- 
tors in the continuous space. In these papers, the authors discussed the critical behavior in 
the thermodynamic limit near the quantum critical point. As a result, the mean- field type 
criticality was confirmed above the upper critical dimension, (e.g., x ~ $ ■) 

III. NUMERICAL VERIFICATION OF MODIFIED FINITE-SIZE SCALING 

In this section, we a ppl y the MFSS to the result of QMC simulation for the single com- 



ponent BH model 18|, [19] . We focus on the superfluid to Mott insulator transition. The 
zero-temperature phase diagram is shown in Fig. [31 which consists of Mott lobes and a su- 
perfluid region. The phase boundary was estimated using the Mott gap. [6( At the tip of the 
Mott lobe, which is a the multicritical point, the dynamical critical exponent z is 1 because 



of the asymptotic particle-hole symmetry. |4j, [l3j The rest of critical lines corresponds to 
the generic transition with the dynamical critical exponent z = 2. In this section, we fix 
the chemical potential as \ij\J = 0.1 and vary the hopping amplitude t/U . Namely, 8 in the 
first argument of the scaling functions corresponds to 5 = t/U — (t/U) c in the present case. 
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We study compressibility k and susceptibility x- Their definitions are 



p 2 dp' 



1 dp 



(5) 



and 



X = - 



1 d 2 F, 
2L d (3 dr] 2 



(6) 



where 



P = 



1 dF 
L d (3 dp ' 



(7) 



The scaling forms of k and \ are derived using the scaling form of the free energy ([T]) as 



We fix the second argument as y = 0.375 and estimate the critical value of t/U as (t/U) c — 
0.088935(7) at p/U = 0.1 using the MFSS of k and x as shown in FigsfTja) and b). In these 
plots, we used the mean-field values for the exponents, leaving the critical value of t/U as 
the only fitting parameter. 

As long as (5t/L 2 = 0(1), we can use the MFSS form just as well as we do in the 
conventional FSS for estimating the critical value of the relevant parameter (t/U in the 
present case). To compare between the estimation using Mott gap and MFSS, we estimate 
the Mott gap at t/U = 0.088935, p/U = 0.5 and plot the corresponding points on the inset 
of Fig. [31 As we see in the figure, the agreement is very good. 

Here, a remark on the range of validity of the MFSS form is appropriate. We consider 
the finite temperature behavior of x a t the quantum critical point 5 = 0. If we neglect the 
applicability condition of the MFSS form and take the limit L — ► oo while keeping f3t finite, 
the finite temperature dependence of x is derived as 



« ~ Y K (x, y) , x ~ L 5/2 Y X (x, y) 



(8) 



where 



x = 5L*l\y 



fit 
L2- 



(9) 




T~ 5/4 (error!), 



(10) 



5 



W r .,W y ,W z ) and 



from the scaling form ([8]). This exponent —5/4 is different from that of mean- field theory 
—3/2. As shown in Sec. \IV\ the reason of this error is that the scaling form (JS]) or (JT]) is 
valid only under the condition of (3t/L 2 = 0(1). To confirm the mean-field exponent, we 
show the finite temperature dependence of x a t the quantum critical point in Fig. [2J 

The superfluid density ps is one of the most important quantity characterizing the su- 
perfluidity. However, it is not straightforward to derive the MFSS form of ps because it is 
not directly obtained from the free energy by simple differentiation. The superfluid den- 
sity ps is proportional to the fluctuation of the winding number W = 
defined as ps = (W 2 ) /((3tL) within the framework of QMC simulation. 
IB] we show that ps = x/{(3L d ), for the model ([2]) in the large- N limit under the condition, 
(3t/L 2 > 0(1), d > 2, and (3t > 1. From the MFSS for x, we obtain, 

p s ~ L-lY ps (x,y). (11) 

Although this form is derived only for the exactly solvable model, we believe that this holds 
in general for the mean-field type critical behavior. We apply this MFSS form to the result 
of ps estimated by QMC simulations. As can be seen in Fig. [T]c), the MFSS (jTTJ) describes 
the data well. 
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IV. LARGE- iV LIMIT OF iV-COMPONENT BOSE-HUBBARD MODEL 

In this section, we consider the model (jSJ) that is known to exhibit a mean-field type 
critical phenomena, to see if the MFSS is applicable to such a model. We consider the 
model on the d dimensional hypercubic lattice in the large- A limit and show that the MFSS 
form Eq. ([I]) is exactly applicable to this case. To derive the self-consistent equation of x i n 
the large- A limit, we represent the partition function as a functional integral by making use 
of a coherent state basis at first. Then, we use the Stratonovitch-Hubbard transformation 
and the saddle-point method,which is also called the steepest descent method. Thus, the 
self consistent equation of susceptibility x in large-A limit is derived exactly as 

i U ^ 1 
X = ~^~t + T7,; r 1 ■ 12 



k exp 



^X- 1 + ¥Eti(l-cosfc, 
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FIG. 1: (Color online) MFSS plots of the single component BH model where [a/U = 0.1, f5t/L 2 = 
0.375. 5 = t/U — (t/U) c with (t/U) c = 0.088935: a) compressibility, b) susceptibility c) superfluid 
density. 
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FIG. 2: (Color online) Temperature dependence of x a t the quantum critical point estimated by 
MFSS (ji/U = 0.1, {t/U) c = 0.088935). The solid line is A(T/t)~ 3 / 2 where A = 0.89. The data 
points are obtained for L = 24, 32, 48. There is no visible size dependence on this scale. The 
statistical error is smaller than the symbol size. 
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FIG. 3: (Color online) Zero-temperature phase diagram of single component BH model in 3D g]. 
FSS indicates the result of MFSS. 
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See Appendix IA II for details of the derivation. By expanding the summand with respect to 



cxp 



, we obtain 



X 



u ^ 



L d 



' L 

J^exp 



n=l 



2v(3t 



cos 



2irn 



-i d 



(13) 



Below we show that this equation has a solution such that x ~ 0((/3i) (/i+2 ' ) / 4 ). Therefore, 
we assume x ~ 0((/3t)^ d+2 '' 4 ) for x i n the r.h.s. of (|T3|) . Then, as shown in Appendix IA 21 
the approximation formula 

L 



E 

i/=i 



J^exp 



n=l 



2z//3t 



cos 



27m 



(14) 



becomes exact in the limit fit — > oo under the condition that d > 2, (3t/L 2 > 0(1). Using 
the self-consistent Eq. ( TTB"]) and the approximation (|14p . we arrive at a simple equation 



X 1 = —fi — t + Ux/3 1 L d . Its solution can be cast into the form, 



z 




t 



L 2 Z 



with a scaling function 



(15) 



(16) 



x + a/x 2 + 1 

At the critical point (p = -t), we obtain xt ~ (/3t) (d+2)/4 x (f3t/L 2 )~ d ^. To make this 
consistent with x — 0(((3t)^ d+2 ^ A ) assumed at first and the condition (5t/ ' L 2 > 0(1), we must 
demand (3t/L 2 = 0(1). Thus, we have proved that Eq. has a solution x = 0((/5t) (d+2)/4 ) 
that satisfies Eq. f|T5|) . and the MFSS form (jSj) has been derived as a formula that is 
asymptotically exact under the condition of d > 2 and f3t/L 2 = 0(1). 

To see the validity of the form of the scaling function (ITBT) . we demonstrate the MFSS 
plot of susceptibility in d = 3. We solve the self-consistent Eq. (fT3"|) without using Eq. (lT4"j) 
and plot on Figs. H]a) and b). As shown in Fig|4]b), the MFSS form fits well in the region 
pt/L 2 = 0(1). 



V. DISCUSSION AND SUMMARY 



In Sec. II III and IIV} we have demonstrated that the MFSS ([T]) is efficient in locating 
quantum critical points whose dynamical critical exponent is z = 2. It has been shown 

9 



a) 




$t/(L 2 Z) 

FIG. 4: (Color online) MFSS plots of susceptibility at U/(t/Z) = 1. a) the x-dependence 
of scaling function of susceptibility P^ Z ^(x,y) and the solution of self-consistent Eq.(|13p with 
y = (5t/{ZL 2 ) = 1, b) y-dependence of P^ Z (x,y) and the solution of self-consistent Eq. (fT3|) at 
quantum critical point 5 = 0. 

that the MFSS is valid only if the second argument of scaling function (3t/L 2 is 0(1). In 
particular, it is not permitted to infinitize L in the scaling forms (jHJ) and (TTTT) while keeping 
(3 fixed. This explains the apparent contradiction between the MFSS and the mean-field 
critical exponents. It should be remarked here that similar situations appear in classical 
models. Suppose that we try to apply the MFSS to a finite-temperature phase transition 
of a classical system and send the system size in some (not all) of the directions to infinity 
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while keeping the size in other directions fixed. Q, [2l|] Singh and Pathria 21] considered a 
system of size L d ~ d ' x L' d where d is larger than the upper critical dimension and d' is less 
than the lower critical dimension. They analyzed a spin model with 0(n) symmetry in the 



limit of V — ► oo. 



21] Then they derived the scaling form of the susceptibility xo as 



2(d-d') ,, /_, 2(d-d') , 

Xo ~ L — Y* [tL— , (17) 



where t = (T - T c )/T c . Then, xo ~ L 2 ( rf - d ')/(4-d') at tne critical point t = 0. On the other 
hand, if we keep L/L' is finite, the MFSS form is 

Xo-L^foiL/L'), (18) 

with the additional argument L/L'. [ijj] If we ignore the validity condition of the MFSS (ITS]) 
and take the limit L' — > oo, we reach an erroneous conclusion, that is xo ~ L d l 2 at t = 0. 

In summary, the MFSS is applied to the quantum critical phenomena with the dynamical 
critical exponent z = 2. Using the A-component BH model, the MFSS form of the sus- 
ceptibility Eq. ffl5l) is exactly derived in the large- N limit with the applicability condition 
d > 2 and f3t/L 2 = 0(1). We also apply the MFSS to the numerical results obtained by 
QMC simulations. As a result, we see that a position of quantum critical point estimated 
by MFSS is identical to that estimated by the Mott gap within the statistical error. Finally, 
note that the scaling function derived in this paper P£(x,y) is in complete agreement with 
the scaling function of 4 model derived in Ref. 9|. While the scaling function is not justi- 
fied by the renormalization group or scaling theories in contrast to the standard FSS below 
the upper critical dimension, the agreement strongly indicates that the mean-field scaling 
function above the upper critical dimension is universal. 
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APPENDIX A: CALCULATION OF LARGE N LIMIT 



1. Self-consistent equation of x 



Here, we derive the self-consistent equation (TX5T) of the Hamiltonian §2§. The partition 
function is expressed as 



Zn 

Sq 



-(So+Si) 



(Al) 



-^E(^(-r)-^(r)+Vl(r)-^(r)) 
Si = ^dr^X)(^(T).^(r)) 2 , 



by the path- integral representation with tpi(r) being an iV-component complex field. Using 
Stratonovitch-Hubbard transformation, the partition function is written as 



P^(T)Pt/>*(r)D Si (T)e- s ° 



x exp 



(W) 



(A2) 



£ {$»(r) (<9 r #(r)) - (/i + ivfo,(r)) C(t)^(t)} 



(A3) 



where s«(r) is an auxiliary field and the integral with respect to Sj(r) is defined as 



(A4) 



In the large- TV limit, using saddle-point method, the auxiliary field Sj(r) is replaced by s 



(see Ref. 



12j and its references ), which makes the exponents of the partition function 



12 



maximum. Using Fourier transformation, we obtain 



Z N = Ae~— S \Z X (s)] 1 



2t 



-pi — iVU s — — cos k 



(A5) 
(A6) 

(A7) 



<5=1 



where A is a some real number caused by the fluctuation of Sj (r) from s, which does not 
contribute to following discussion and the product of k runs over the first Brillouin zone 
k = (27r/L)(rai, • • • , ma), with rrii = 1, 2, • • • , L. The stationary solution s must satisfy 



d_ 

ds 







0. 



(A8) 



which yields, 



L d 



(A9) 



k 



The susceptibility \ is related to s by 

1 rfi 



X 



X 



dr ^>*(r) ■ MO)) = (-fi-t- iVUs 



(A10) 



Therefore, \ satisfies 



U 



X 



u s. — ■> 



fc exp 



¥ Etia-cos^+^x- 1 



(All) 



2. Derivation of Eq.flHD 



In Sec. IIV| we derive x — O ((/3i/Z)( 2+d )/ 4 ) by self-consistent analysis. Namely, assuming 
the condition x = O ((/%/Z)( 2+d )/ 4 ), we prove the resulting solution satisfying this condition. 
Here, assuming 



(f3t/Z)M/* x -i = 0(1), 



(A12) 



we provide an approximation form 

L 



J^exp 

n=l 



~z- r- cos hr 



(A13) 
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which becomes exact under the condition that 



d > 2, 
(3t/L 2 > 0(1), 



(A14) 
(A15) 



and 



Pt > 1. 



(A16) 



To begin with, we rewrite the l.h.s. as 



u=l 



^exp 



n=l 



2vBt f /2vm 
1 — cos 



S>2 e"^ 1 [1 + A 



d 



i/=i 



i/=i 



^ a\(d-a)\ v 

a=l v ' 



(A17) 



where 



-l+L/2 

A v = e^ l ' z + 2 ex P 



n=l 



2vQt |\ /2vm 
I — cos I — — 



(A18) 



Here we note that (3x~ l — 0. (This is because /?x 1 = 0((/3i)~' d ~ 2 -'' 4 ) (by the condition 
Eq. (lA12p ) and by the condition Eq.( lA14l) this is vanishing in the limit of Eq.( 1A16l) .) Since 
f3\ l — 0, the first term of the r.h.s. of Eq. (lA17p is approximated by the formula 



oo 

^e"^" 1 = p- x X + ((/3X- 1 ) ) = O {{fit) 



d-2 
4 



(A19) 



u=l 



Below we show that the second term of Eq. (1A17j) is a correction term that vanishes in the 
limit of f3t — > oo. At first, A v is bounded as 



L/2 

< A v < 2 exp 

n=l 
rL/2 

< 2 / dp exp 

Jo 

rL/2 

< 2 / dp exp 

Jo 



2v/3t f /2im 
1 — cos 



2u/3t L /2vrp 
1 — cos 



2upt f8p 2 



< 



nL 2 Z 
16V/3? 



(A20) 
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Then, the second term of Eq. (1A17I) is evaluated as 
d 



d\ 



° < S a !(j-n)! 

a=l v ' 



< 



d! 



n! -n)! 

tZ! 



,f=i 

oo 

E 



" ^-f a! (d-a)! 1 16/% 



f7TL 2 Z 

Vovfit 



a/2 



a/2 r 



E 



d! 



^a!(d-a)! [l6(/5t/L 



7Tz 



a/2 



dge 9 g 



(ftT 1 ) 



ot-2 
2 



(A21) 



Since j3t/L 2 > 0(1), and /?x 1 ^ 1) the a = 1 term is dominant. Therefore, the sec- 
ond term is of the same order as ((3t/ X 2 ) -1 / 2 ^- 1 )- 1 ^ By the condition Eq.flSG), this 
is O ((/3t/L 2 ) -1 / 2 x (/?t)( d ~ 2 ^ 8 ). Therefore, the ratio of the second and the first term of 
Eq. (jA17p becomes less than O ((fit/L 2 ) -1 / 2 x (f3t)~( d ~ 2 ^ 8 ) . This is vanishing because of the 
condition Eqs. (TAHI) . fTATBD and (TATdI) . Thus Eq. flAT3D has been derived. 



APPENDIX B: SCALING FUNCTION OF SUPERFLUID DENSITY IN LARGE- 
N LIMIT 

In this section, we provide that the MFSS form of superfluid density using the N- 
component BH model. The outline of this section is as follows. First, we obtain the explicit 
definition of superfluid density, which is estimated using the winding number in QMC sim- 
ulation, with an infinitesimal twist of phase of bosonic operator. Next, we calculate the 
superfluid density of A^-component BH model exactly. The result reveals that the superfluid 
density ps is proportional to the susceptibility x- Then, we derive the MFSS form of ps as 
that of X- 

To start with, we derive an expression for the superfluid density ps introducing an in- 
finitesimal twist of phase of bosonic operators. Namely, we modify the Hamiltonian (J2J) by 
b a i ~ *■ ^ a i e%er ^ i bai b ai e~ l0r i , where rf is the ^-coordinate of the site i. (Because of the 
periodic boundary condition, 9 should be discrete. That is, 9 = 2irn/L where n is integer. 
However, considering a sufficiently large system, we regard 9 as a continuous real number.) 
Then, we define the twisted Hamiltonian H.N6, the partition function Z^g and the free energy 
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Fno as, 



N N 

a=l a=l i 

N N 

+^EEE 6 - 6 kM«* (bi) 

a=l (3=1 i 

Z m = Tr [e'^o] , (B2) 
-FVe = ~\nZ m . (B3) 



The superfluid density p s is denned with this twisted free energy F Ne as, 

1 d 2 F\ 



NO 

PS ~ 



(B4) 



0=0 



2f3L d (t/Z) d6 2 

Next, we calculate the ps in the large- N limit. The partition function Z^e is obtained as 
well as the non-twisted partition function (See Appendix. IA II) as, 

Z m = Ae-^ 2 [Z w (s)] N , (B5) 
fc 

2t t\ 2t 

Afcg = —p — iy/Us — —\^cos ks — — cos (k z + 6) . (B7) 

Zj Z 

(5=1 

The derivation of the free energy is straightforward using this partition function. Then, we 
obtain the superfluid density is 

k 

with Afc defined in Eq. flA7j) . This superfluid density is smaller than the total density of 
particle, 

1 8Fnq = q 



P = 



L d 

k 

and larger than the density of particles of k = 0, 

1 1 



L d p dp ' 

y E e /3A fe _ x ' ( B9 ) 



* = WTSF^T- (B10) 
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That is, 

Po < ps < P- 

As shown in Appendix IA 2\ 



X 

P = Po 



L d f3' 

under the condition (3t/L 2 > 0(1), d > 2 and (3t — > oo. Using the inequality (IBllj) 



obtain 

X 



Ps 



L d p' 

As shown in Sec. IIVt we derive the MFSS form of ps as 



PI (x, y) = 



x + \j x 1 + Auy 1 

The applicability condition of this MFSS form is d > 2 and f3t/L 2 = 0(1). 
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